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ABSTRACT 

Aims. We investigate the formation of neutral and singly ionized scandium lines in the solar photospheres. The research is aimed 
derive solar logg/£ Q (Sc) values for scandium lines, which will later be used in differential abundance analyses of metal-poor stars. 
Methods. Extensive statistical equilibrium calculations were carried out for a model atom, which comprises 92 terms for Sc i and 79 
for Sc n. Photoionization cross-sections are assumed to be hydrogenic. Synthetic line profiles calculated from the level populations 
according to the NLTE departure coefficients were compared with the observed solar spectral atlas. Hyperfine structure (HFS) broad- 
ening is taken into account. 

Results. The statistical equilibrium of scandium is dominated by a strong underpopulation of Sc i caused by missing strong lines. It 
is nearly unaffected by the variation in interaction parameters and only marginally sensitive to the choice of the solar atmospheric 
model. Abundance determinations using the ODF model lead to a solar Sc abundance of between log£ Q = 3.07 and 3.13, depending 
on the choice of / values. The long known difference between photospheric and meteoritic scandium abundances is confirmed for the 
experimental /-values. 
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1. Introduction 

On ! 

CD ■ According to the theory of nucleosynthesis, the a elements 
' are mostly produced by Type II supernovae, while some iron- 
, peak elements have significant contributions from Type la su- 
£NJ ■ pernovae. The synthesis process and sites of scandium, as an el- 
CD ement intermediate between a elements and iron-peak elements 
OO , in the periodic table, are not clear now. 

The variation of the scandium abundance pattern in long- 
^ , lived F- and G-type stars with different metallicity can provide 
k> ' some information on the element nucleosynthesis and the chem- 
rN \ ical evolution of the Galaxy. There is an unresolved inconsis- 
■ tency between different Sc abundance results. In some analyses, 
an Sc enhancement relative to Fe is found in metal-poor stars 
(e.g. Zhao & Magain 1990 Nissen et al. 120001 1; however, others 
(e.g. Gratton & Sneden [19911 Prochaska & Mc William [20001 
found no evidence of any deviation from [Sc/Fe] = 0.0. 

Generally, the solar photospheric abundances serve as a ref- 
erence for abundance determinations in metal-poor stars, so a 
reliable set of photospheric abundances is important. Ever since 
Anders & Grevesse ( 1989 ) published their widely used solar el- 
emental abundance table, many revisions and updates to pho- 
tospheric and meteoritic abundances of the elements have be- 
come available, although the solar photospheric scandium abun- 
dance has not been updated for quite a long time. The photo- 
spheric abundance value ofloge (Sc) = 3.10 + 0.09 adopted by 
Grevesse ( 1984) was changed to 3.05+0.08 by Youssef & Amer 
( [19891 . Neuforge ([19931 obtained 3.14 + 0.12 from the Sc i lines 
and 3.20 + 0.07 from Sc n lines. The average value of 3.17 +0.10 
was adopted by Grevesse & Noels d 1993b and was kept in the 
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newest tabular version of Grevesse et al. d2007l l. which is some- 
what higher than the meteoritic value of 3.04 ± 0.04. 

It should be noted that local thermodynamic equilibrium 
(LTE) has been assumed in previous papers about scandium 
abundance determinations, and NLTE investigation of the scan- 
dium element has never been published. In general, departures 
from LTE are commonplace and often quite important, in partic- 
ular for low surface gravities or metallicities, with minority ions 
and low-excitation transitions the most vulnerable (see the re- 
view paper of Asplund[2005). In the Sun and other near-turnoff 
stars, the ionization energy of Sc i (6.56 eV, see Fig.[TJ indicates 
that this is a minority ion. Lines of Sc i should therefore be more 
susceptible to NLTE, because any small change in the ionization 
rates largely changes the populations of the minor ion, although 
there is no guarantee that Sc n lines behave properly in a standard 
LTE analysis. 

In this paper we investigate the statistical equilibrium and 
formation of neutral and singly-ionized scandium lines in the 
solar photosphere. We note that we do not intend to study the 
influence of atmospheric inhomogeneities on any scale or that 
of chromospheric temperatures and pressures. The method of 
NLTE calculations is briefly introduced in Sect. 2. In Sect. 3 the 
synthesis of the Sc i and Sc n lines under NLTE and hyperfine 
structure is presented. The discussion fills the last section. 

2. NLTE line formation calculations 

Abandoning the LTE approximation introduces a great deal of 
additional complexity into the line formation calculations. Under 
NLTE conditions, the atomic populations are described by a set 
of statistical equilibrium equations in which radiative and colli- 
sional processes are to be taken into account. Our calculations 
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Fig. 1. Adopted level structure of the Sci and Sen atom taken from the NIST data bank. Explicitly calculated line transitions are 
indicated, where continuous lines refer to allowed and dotted lines to forbidden transitions according to selection rules assuming LS 
coupling. 



were carried out with a revised version of the DETAIL program 
(Butler & Giddings 119851 ), which solves the radiative transfer 
and statistical equilibrium equations by the method of acceler- 
ated lambda iteration. 

Our calculations were performed on a partially ionized back- 
ground medium consisting of a plane-parallel, homogeneous, 
line-blanketed theoretical model of the solar photosphere, which 
includes a simple approach to convective equilibrium based on 
the mixing-length theory. The model was computed with the 
MAFAGS code (Fuhrmann et al. 11997ft . In contrast to the line 
formation itself, the atmospheric model assumes LTE to obtain 
the final temperature-pressure stratification. It uses opacity dis- 
tribution functions (ODF) for line-blanketing, based on Kurucz 
(1992), and calculated with opacities rescaled to a solar iron 
abundance log e Q (Fe) = 7.51 (more details are found in Gehren 
et al. |2001| ). The resulting atmospheric stratifications of tempera- 



ture and pressure are similar to those given by other solar models 
(see comparison in Grupp 2004 Fig. 15). 

2.1. Atomic model 

A comprehensive atomic model is required for NLTE calcula- 
tions . Similar to other iron-group elements, scandium has a 
complex atomic structure. Our atomic reference model is con- 
structed from 256 and 148 levels for neutral and singly-ionized 
scandium, respectively. Energies for these levels are taken from 
the NIST data banlfl. After a few early test calculations with this 
complete fine structure model, we found that the calculations 
could be considerably reduced by combining all fine structure 
levels into 92 and 79 terms for Sc i and Sc n, respectively. The 
corresponding fine structure data were appropriately weighted 
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to determine term energies. The atomic term model used for our 
calculations is displayed in Fig.Q] It shows that completeness is 
fading at high excitation energies, with energy gaps of between 
0.3 and 1.0 eV for the neutral doublets and quartets, and gaps of 
~ 1.5 eV for the ionized singlets and triplets. 

The number of bound-bound transitions treated in the 
NLTE calculations is 1104 for Sci and 1034 for Sen, numbers 
again considerably reduced from the original level transitions. 
Wavelengths and oscillator strengths of the fine structure transi- 
tions are taken from Kurucz's database (Kurucz & Bell 1995 ), 
and they are weighted by statistical weights to yield artificial 
term transitions. A Grotrian diagram for Sc i and Sc n is dis- 
played in Fig. Q] Solid and dotted lines represent the allowed 
and forbidden transitions included in the model atom, respec- 
tively. Since the reduction of the fine structure model atom does 
not affect the resulting calculations of the population densities, 
it is unlikely that the hyperfine structure (HFS) has any direct in- 
fluence on the NLTE results, particularly since the known HFS 
is small. 

For bound-free radiative transitions in the Sc atom, 
hydrogen-like photoionization cross-sections (Mihalas |1978t are 
adopted, because data from the Opacity Project (OP; see Seaton 
et al. 119941 ) are not available. In our current analysis this may 
be the most uncertain representation. In previous studies of Fe 
(Gehren et al. 2001 ), and K (Zhang et al. 120061) . where complex 
calculations of such cross-sections were available, we found that 
hydrogenic approximations occasionally tend to underestimate 
the photoionization cross-sections by one or two orders of mag- 
nitude. The effect on the NLTE analysis is examined below. 

As usual, background opacities are calculated with an opac- 
ity sampling code based on the line lists made available by 
Kurucz ( 1992). Since background opacities affect the photoion- 
ization rates directly, their consideration is important. We note, 
however, that the millions of faint lines, which may be somewhat 
more important for model atmosphere construction, are marginal 
for our line formation calculations. 

In our calculations for Sc, we take into account inelas- 
tic collisions with electrons and hydrogen atoms leading to 
both excitation and ionization. Because laboratory measure- 
ments and detailed quantum mechanical calculations for col- 
lision cross-sections are absent, approximate formulae are ap- 
plied. The formulae of van Regemorter ( 1 19621 ) and Allen ( 119731 ) 
are used to describe the excitation of allowed and forbidden 
transitions by electron collisions, respectively. Ionization cross- 
sections for electron collisions are calculated with the formula 
of Seaton d 19621 ). Drawin's d 19681 119691 ) formula as described 
by Steenbock & Holweger (1984) is used to calculate neutral 
hydrogen collisions. Recently it was indicated both experimen- 
tally and theoretically that Drawin's formula overestimates the H 
collisional cross-section by one to six orders of magnitude (e.g. 
Belyaev et al. 1999 Belyaev & Barklem 20031. so a scaling fac- 
tor Sh is applied to the Drawin formula in our calculations, for 
which results are given below. 

2.2. NLTE level populations 

The atomic model described above still has a number of free 
parameters that represent our ignorance of the microscopic in- 
teraction processes. Whereas the number of levels (terms) and 
lines (transitions) comprises the basic structure of the two lower 
scandium ions, it is the interactions that require some additional 
fine tuning. As explained above, we introduce standard multi- 
plication factors for electron collisions, hydrogen collisions, and 
photoionization cross-sections. These factors are defined with 
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Fig. 2. Departure coefficients b„ for some levels of Sci and 
Sc ii in the ODF solar atmospheric model applying our standard 
atomic model with multiplication factors S & - 1, Sh= 0.1, and Sp= 
1 for bound-bound electron collisions, bound-bound hydrogen 
collisions, and bound-free photoionization coefficients, respec- 
tively. The kinetic equilibrium calculations are described in the 
text. 



respect to the standard formulae for these three types of pro- 
cesses, which were mentioned in section 12.11 Using the stan- 
dard atomic model, departure coefficients fcj = «^ LTE /n^ TE for 
Sc i and Sc n terms in the solar atmosphere are presented in Fig. 
|2] The standard modelrefers to 5 e = 1, Sh- 0.1, and Sp= 1. 

Since there are no other strong indicators, it is necessary to 
explore the influence of such a parametric variation in interac- 
tion strengths on the level population densities. This is done by 
varying only one parameter at a time while holding all others 
at their standard values, and the results are surprising in that 
we find no strong influence for either of the multiplication fac- 
tors (Fig. [3]). All calculations represented here document the ex- 
treme weakness of the solar Sc i lines resulting from both the 
low abundance and the low ionization potential. All lines thus 
form in the same atmospheric depths as the local FT contin- 
uum. There the mean integrated line intensities become con- 
stant (and the net bound-bound radiative rates positive, or up- 
wards) outside r =; 1 . All the lower levels (those with less than 
3 . . . 4 eV excitation energy) are therefore depopulated by this 
collective pumping process, whereas the upper levels are simul- 
taneously populated. However, no true population inversion is 
achieved. Photoionization rates from the lower levels are low, but 
mostly higher than the corresponding recombination rates. This 
again supports the depopulation of the lower levels throughout 
the solar atmosphere. Of course, the radiative rates are modified 
by collision rates, and the total rates are driven more towards 
zero net rates, but this does not prevent the general depopulation 
trend. 

Figure|3]demonstrates clearly that even relatively large vari- 
ations in the multiplication factors do not change the run of the 
level populations too much. There are thermalizing effects when 
increasing electron or hydrogen collision factors, but these affect 
the population ratios, and not the departure coefficients them- 
selves very much. A notable exception is represented by the S e = 
10, because this starts to couple the lower terms more efficiently 
to the upper ones and thus thermalizes the whole Sc i system. 
The lack of variation with Sp is simply due to the decoupling of 
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Fig. 3. Departure coefficients b n for levels of Sc i and Sc n in the ODF solar atmospheric model testing the multiplication factors S e = 
0.1 and 10 (top), Sh= and 1 (middle), and 5p= 0.1 and 10 (bottom). 



the radiation field and the small hydrogen-like photoionization 
cross-sections. 

Sc n is the dominant ion of the element under solar atmo- 
spheric conditions, with more than 99% of the scandium atoms 
being ionized under the atmospheric conditions found in the Sun. 
Its lines are substantially stronger than those of Sc i, although not 
comparable in line strength with other metals. Only the strongest 
lines form outside logr -1. Figure [3] shows that all Sen lines 



in the visible spectrum are formed near to LTE conditions. There 
is no clear indication how to select the proper scaling factors. In 
view of the minor changes due to parameter variation, we choose 
typical factors to establish the standard model atom (see Fig. [2j. 
Since this choice may affect the line formation in stars different 
from the Sun, we will extend the analysis for metal-poor stars to 
obtain a more significant choice of the scaling factors. 
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Table 1. Atomic data and abundance results for Sc i and Sc n lines. 
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3. Analysis of scandium lines in the solar spectrum 

In this section, we investigate the formation of Sc i and Sc n 
lines in the solar atmosphere and derive the scandium abundance 
based on spectrum synthesis. Lines in the solar spectrum are cal- 
culated using the plane-parallel hydrostatic MAFAGS-ODF so- 
lar model atmosphere with r eff = 5780 K, logg = 4.44, [Fe/H] 
= 0.00, £ t = 0.90 km s _1 (for a more detailed comparison with 
opacity sampling models see Grupp l2004l l. An initial scandium 
abundance of log e (Sc) = 2.99 is adopted here. It should be 
noted that the atmospheric model of the Sun does not depend on 
the scandium abundance. For all elements except scandium, we 
assume LTE. 



0.15 r 



-0.05 |_ j 
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log gf (Kurucz) 

Fig. 4. Comparison of logg/-values from different sources. 
Filled circles represent the Sc i ion, open circles refer to Sc n. 
A log gf means Kurucz - Lawler & Dakin. 



3.1. Atomic line data 

For the solar abundance analysis we selected 4 Sc i and 17 Sc n 
lines, which ideally should satisfy the following conditions: they 
are relatively free from blends, and oscillator strengths and hy- 
perfine splitting parameters are available. However, both condi- 
tions are not always guaranteed. In particular, HFS data for the 
ground state of Sc n are missing making the analysis of the cor- 
responding lines more uncertain. Unfortunately, the number of 
sufficiently strong lines in the visible of both ions is very lim- 
ited, and in metal-poor stars this requires concentration on the 
leading lines of Sc n, since the equivalent width of Sc i 4023 in 
typical turnoff stars is less than 0.5 mA. 

To determine the solar abundance, it is necessary to know 
the accurate oscillator strengths (/ values) of the spectral lines. 
Two sets of oscillator strengths are applied and compared in 
our abundance determinations: (i) theoretical values taken from 
Kurucz' databasfl and (ii) experimental data of Lawler & Dakin 

http://kurucz.harvard.edu, see also Kurucz & Bell JT9951 



( 119891 ), which were based on lifetimes measured by Marsden 
et al. ( 119881 1 together with branching fractions. Van der Waals 
damping constants log are computed according to the Anstee 
& O'Mara ( 119911 1 19951 1 interpolation tables. Input parameters 
needed to perform spectrum synthesis for the selected lines are 
given in Table 1. The different sets of /-values from Kurucz 
and from Lawler & Dakin are compared in Fig. [4] There is a 
small systematic offset between the two data sets, with a mean 
A log gf (Km - LD) 0.05. The corresponding line data for our 
line fits are reproduced in Table 1. 

In solar system matter, scandium is represented only by 
the 45 Sc isotope. Similar to other odd-Z elements, hyperfine 
structure interactions between nuclear and electronic wave func- 
tions split the absorption lines of Sc into multiple components. 
Hyperfine structure components of line transitions are calcu- 
lated as usual from magnetic dipole splitting constants, A(J), 
and electric quadrupole splitting constants, B(7), of the corre- 
sponding levels. For most of the lines the HFS components fall 
into small intervals; we therefore combine all components within 
5 mA (Sc i) or 10 mA (Sc n), which reduces many HFS patterns 
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Table 2. Hyperfine structure data for Sc i and Sc n levels used in 
this analysis. 
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0.0209 0.00363 -0.00124 Childs (1971) 

1.4395 0.0083 -0.0003 Ertmer & Hofer (1976) 

1.4478 0.0095 -0.0004 Ertmer & Hofer (1976) 

3.6191 0.0029 Ertmer & Hofer (1976) 

3.6332 0.0015 Ertmer & Hofer (1976) 
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in units of cm 



to two or three coadded lines. The basic data are given in Table 
[2] Abundance differences with respect to the full HFS pattern are 
all within 0.01 dex. 



3.2. Line profile fitting 

The observed solar flux spectrum was taken from the Kitt Peak 
Atlas (Kurucz et al. 1984). Spectrum synthesis was employed to 
determine the abundance of scandium in the solar atmosphere. 
As in earlier work we used the interactive spectral line-profile 
fitting program SIU, which is an IDL/Fortran software package 
(Reetz |T9"9"0"l l. To match the observed spectral lines, the synthetic 
spectra were convolved with a mean solar rotational velocity of 
1.8 km s _1 and a radial-tangential macroturbulence S rt which 
is found to vary for lines of different mean depth of formation 
between 2.8 and 4.0 km s . 

Except for the obvious influence of strong metal or Balmer 
line wings in the spectral range under consideration, we reset 
the local continuum position interactively to the maximum flux 
in a +5 A interval around the line center. This is never a problem, 
because our spectrum synthesis includes all important lines and 
thus allows confirmation of the local maximum flux estimate. 
Any uncertainty in this process, even in the yellow wavelength 
range between 5700 and 5850 A where the solar atlas displays 
a continuum that is systematically high by ~ 2%, is smaller than 
the general profile fitting error. Our estimate of its influence on 
the abundances is ~ 0.01 dex. 

Line profiles are computed under both LTE and NLTE as- 
sumptions: fitted to the observed profiles by means of scandium 
abundance variations. Column 10 of Table[T]reproduces the log- 



arithmic abundances log e of the fits including a number of weak 
blends on either side of the profiles. They are based on Kurucz' 
log gf values and NLTE level populations. The logarithmic cor- 
rections due to weak blends are listed in the Ab column. Thus, 
fitting the lines without blends would have resulted in a higher 
logarithmic abundance log s - Ab. The difference required to fit 
LTE and NLTE profiles is referred to as the NLTE correction 
( AX = loge NLTE - loge LTE ); it is given in column 13 of Table Q] 
for each line. Equivalent widths from NLTE profile fits (see Fig. 
|5]l are given in column 6 of Table Q] The last three columns give 
the number of components, the maximum wavelength separation 
(in mA) of the hyperfine structure (HFS) lines and an asterisk, if 
HFS data are missing for one of the levels(cf. Table|2]). 

Some of the synthetic profiles for selected lines, together 
with the observed solar spectrum, are presented in Fig. [5] 
Whenever available, we have included known blends. For com- 
parison we show profiles under LTE and NLTE conditions in- 
cluding HFS. Generally, all neutral lines are much fainter in 
NLTE due to significant underpopulation, whereas lines of ion- 
ized scandium are slightly stronger under NLTE conditions. 

A marginal line core asymmetry with the observed line bisec- 
tors shifted 1 or 2 mA to the red seems to exist. This is known 
from the solar spectrum synthesis of other metals, such as Si, 
Ca, or Mn. It is probably the result of hydrodynamic stream- 
ing patterns that cannot be represented by our simple micro- 
/macroturbulence scheme. A more disturbing defect is the near- 
continuum flux deficiency of the red line wing, best seen in 
the 5657 and 6245 A lines. This feature is also found in other 
metals. Sometimes addressed as the result of a weak line haze 
or blends, the systematic blue-red asymmetry is more likely to 
result from hydrodynamic flows, too. We emphasize, however, 
that this wavelength region displays a particularly disturbed so- 
lar continuum flux. Table Q] already gives a hint that some 
of the lines synthesized here suffer from missing HFS data. 
Unfortunately, the strongest lines of both ions are affected, re- 
quiring more detailed comments. 

Sci 4023.69 A has been calculated with only the HFS split 
of the a 2 D 5 ^ 2 level. That this is possibly a fair but not perfect 
representation of the true hyperfine structure width of the line 
is documented in Fig. [5^, where the synthesized line halfwidth 
fits that of the solar spectrum. However, the total separation of 
the HFS components is small (see TableQ}. Another distortion is 
caused by a number of faint line blends on the red core and wing 
of this line. The known components are two highly excited faint 
Mm lines at 4023.72 and 4023.84 A, for which no HFS data 
are available. There is also a faint Cri line at 4023.74 A. These 
blends have been considered in Fig. [5^, but even the two lines 
within 50 mA of the Sc i line center have no influence on the fit 
of the core. Introducing these blend components and fitting the 
full core profile results in a Sc line abundance change below 0.01 
dex. The remaining red wing depression centered on 4023.84 A 
may be the result of the unknown Mn i HFS. 

Sci 5081.55, 5671.80, and 5686.83 A are substantially 
fainter neutral Sc lines, with solar equivalent widths of only 12, 
14, and 9 mA, respectively. That makes the line abundances be 
sensitive to the continuum position, which is particularly uncer- 
tain between 5600 and 5700 A. Only one of the lines is presented 
in Fig. [5J). Although most of the Sc lines were chosen to be as 
free of blends as possible, the wings of the 5671.80 A line are 
covered by quite a few weak lines of Mn i, Fe i, and Ti i, all of 
which are in the range of equivalent widths between 0.5 and 1 .5 
mA. These lines have central depths between 1 and 2% of the 
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Fig. 5. Representative solar line profiles for Sc i and Sc n. Always shown are the NLTE and LTE profiles for the given abundance, 
where the NLTE profiles refer to the best fit. A 4246.83 additionally shows the NLTE line with the originally calculated HFS 
components (second component at 4246.839 instead of 4246.844 A, see text). 



continuum flux, and their (gaussian) wings give some combined 
contribution to the Sc line core. Thus the full influence of the 
weak blends on the Sc abundance of this line is -0.04 dex. For 
the other two lines, the blend corrections are -0.04 and -0.07, 
respectively. 

Sen 4246.83 A is one of the strongest lines of Sen in the 
visible. It contains two weak blends on the blue line wing. Due to 
missing HFS data for the lower level, a l D 2 , the resulting profiles 



are not as reliable as those obtained for the excited levels. This is 
evident from Fig. [5};, where the line width of the profile with the 
originally calculated HFS components does not fit the observed 
solar profile unless the HFS separation is increased by 5 mA 
moving the second component from 4246.839 A to 4246.844 A. 
This is a purely empirical correction; however, it is the only way 
to fit the solar line profile without postulating an unknown blend. 
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Fig. 6. Solar Sc abundances based on Kurucz gf values for lines 
of Sc i (filled circles) and Sc n (open circles) plotted as func- 
tion of the line equivalent widths. The mean abundance from 
Sc i lines is represented by the dotted line, that of Sc n by a con- 
tinuous line. Top: LTE, bottom: NLTE. 



Sen 5526.81 A and 5657.87 A (Figs.EJi and e) represent the 
lines on the flat part of the curve-of-growth; i.e., they strongly 
depend on the microturbulence parameter. Their hyperfine struc- 
ture is known and seems to fit the solar spectrum for both lines. 
There are two blends on the red wing of the 5526 A line, but not 
near enough to the core of the Sc line to affect the abundance. A 
few weak neutral Cr, Fe, and V lines on both wings of the 5657 
A line have no influence on the abundance either. 

Sen 6245.63 A belongs to the group of well-separated, un- 
blended weak lines. Again, the hyperfine structure fits the ob- 
served spectra as shown in Fig.|5f. This line lies on the extended 
wing of Fei (816) 6246.334 A, but that does not change the line 
abundance. 

Altogether, the remaining influence of blends is small, with 
Ab = -0.021, but it is systematic in that it always reduces the 
mean Sc abundance. It is even more important because it tends 
to reduce the peaks of the abundance distribution. 

3.3. The solar scandium abundance 

Our method of spectrum synthesis yields the product of the os- 
cillator strength for a given transition and the abundance of the 
element, log gfs e . The results are reproduced in Table 1 assum- 



ing NLTE conditions. Using the values obtained for log gfe Q 
and the log gf values from different data sets, we computed Sc 
abundances for the individual lines. Figure 6 shows LTE and 
NLTE abundance results based on the Kurucz calculated oscilla- 
tor strengths for all lines as a function of their equivalent widths. 

Under the LTE assumption, absolute solar abundances de- 
termined from Sc i lines are significantly lower than the values 
obtained for the Sc n lines. The mean LTE abundances for 4 Sc i 
and 17 Sen lines are 2.90 + 0.09 dex and 3.10 + 0.05 dex, respec- 
tively. This discrepancy of the ionization equilibrium is resolved 
in NLTE calculations. Under NLTE assumption, Sc i and Sc n 
lines give very consistent abundance results, i.e. 3.08 + 0.05 dex 
and 3.07 + 0.04 dex, respectively. Using Kurucz' gf values, the 
mean abundance for all 21 Sc lines under NLTE is 3.07 ± 0.04 
dex. 

Using instead the laboratory log gf values of Lawler & 
Dakin, the mean abundance of 17 Sc lines under NLTE is 
3.13 + 0.05 dex. The ionization equilibrium differs by 0.03 dex. 
While the results for both sets of gf values are a brilliant justi- 
fication of the NLTE assumption itself, it seems that the experi- 
mental gf values provided by Lawler & Dakin ( 1989) should be 
slightly more reliable than those calculated by Kurucz. However, 
the standard deviation of the experimental data is slightly higher. 

4. Discussion 

The most interesting result of this analysis turns out to be that 
scandium is the first element found to show strong solar NLTE 
abundance effects in one of its ions. This appears to be the con- 
sequence of missing strong lines in Sc i. All the lines of this ion 
are weak because of two contributions: (a) the solar scandium 
abundance is 3 dex lower than that of other metals, and (b) due 
to its low ionization energy of only 6.6 eV, Sc i is an extreme mi- 
nority ion, that is 99.8% ionized in its solar region of line forma- 
tion. Therefore it does not help that the gf values are relatively 
normal, and the statistical equilibrium of Sc i is far from ther- 
mal even in layers of continuum formation. The departure co- 
efficients calculated for different values of photoionization and 
collision parameters (see Fig. |3) reveal the insensitivity of our 
results to the choice of details of the model atom. This may no 
longer be the case in more metal-poor stars. 

Although still far from complete, the experimental hyperfine 
structure data are compatible with the observed solar line pro- 
files, with the exception of a few lines, for which either lower or 
upper level HFS data are missing. We have not attempted to de- 
termine abundances without HFS, because the profile fits were 
far from realistic in most cases (in particular for Sc n). The gen- 
eral trend of such an analysis would be a relative increase in 
the abundances as a response to a decrease in line broadening, 
which affects mostly the few lines on the flat part of the curve- 
of-growth. 

Depending on whether calculated gf values of Kurucz & 
Bell d 19951 1 or experimental gf values of Lawler & Dakin (1989) 
are preferred, the solar scandium abundance is 



log e Sc , = 



3.07 + 0.04 
3.13+0.05 



(Kurucz & Bell) 
(Lawler & Dakin) . 



For the experimental gf values, the difference with respect 
to the meteoritic scandium abundance turns out to be as high as 
before. Therefore the emphasis lies on our current analysis be- 
ing much more restrictive than recent results. If the higher pho- 
tospheric abundance were caused by a systematic error in our 
analysis, this should be found in Fig. [6p. Assuming say that all 
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lines were affected by remaining unknown blends, the resulting 
abundances should show a trend with line strength (weak line 
abundances are more affected than strong lines). Since such a 
trend is not found, we may rule out the importance of unknown 
blends. A similar argument holds for the role of the HFS, how- 
ever, now producing an inverse trend, where weak lines should 
not depend on the HFS. 

A completely different systematic abundance error could re- 
sult from the solar model atmosphere. In the above investigation 
we started the statistical equilibrium and the synthetic spectrum 
calculations with our standard ODF model atmosphereQ, which 
is very much the same as that of Kurucz ( 1992), and it makes use 
of his opacity distribution functions, corrected for metal abun- 
dance. Such ODF models tend to have slightly lower temper- 
atures near optical depth t =s 1. Opacity sampling (OS) mod- 
els are different. Their higher temperatures are essentially the 
reason for the insufficient fit the solar Balmer lines (see Grupp 
2004). Since the temperature difference for the solar OS model 
of Grupp is only around 40 K or less, it changes the Sc line for- 
mation by a negligible amount. The abundance entries in TableQ] 
document that there is a mean abundance difference in the sense 
ODF-OS of only -0.016 dex. 

Comparison with the results obtained for the ODF model at- 
mosphere shows that the solar scandium abundance for the lab- 
oratory gf values of Lawler & Dakin, logesc.o = 3.13 + 0.05, 
is off the meteoritic value of 3.04 by 0.09 dex. A large fraction 
of the remaining scatter of the single line abundances seen in 
Fig. |6j) is probably caused by the uncertain gf values, both for 
calculated and experimental data. While previous analyses could 
not confirm a correspondence between photospheric and mete- 
oritic Sc abundances due to a relatively large line-by-line scatter, 
our results reduce the problem to a simple discrepancy. Kurucz' 
/-values lead to a solar Sc abundance well in agreement with 
the meteoritic value, whereas the experimental data of Lawler 
& Dakin deviate from that reference by nearly 2<x. Currently, 
there is no hint as to why photospheric scandium should differ 
in abundance from that found in chondrites. 
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